Pair correlation function for spin glasses 
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We extract a pair correlation function (PCF) from probability distributions of the spin-overlap 
parameter q. The distributions come from Monte Carlo simulations. A measure, w, of the thermal 
fluctuations of magnetic patterns follows from the PCFs. We also obtain rms deviations (over 
different system samples) Sp away from average probabilities for q. For the linear system-sizes L we 
have studied, (i) w and 5p are independent of L in the Edwards- Anderson model but scale as l/L 
and y/L, respectively, in the Sherrington-Kirkpatrik model. 

PACS numbers: 75.10.Nr, 75.30.Kz, 89.75.-k, 75.40.Mg 



Introduction — After decades of much work, the nature 
of the spin-glass (SG) phase is still unclear. A SG phase 
is found at low temperature in magnetic systems which 
are both quench- disordered and frustrated Random 
spin positions as well as random spin-spin couplings are 
sources of quenched disorder. When competition arises, 
because not all spin-spin coupling energies can be mini- 
mized simultaneously, a system is said to be frustrated^ 
Fixed disorder and built in competition are the two es- 
sential ingredients of complex systems 

The Sherrington-Kirkpatrick (SK) spin-glass modelj^ 
in which each spin-spin coupling is assigned at ran- 
dom, without regard to spin-spin distance, is quench- 
disordered and frustrated. Its exact solutioni*^"— im- 
plies different random number seeds (which fully spec- 
ify all couplings) can give rise, in thermal equilibrium, 
to magnetic patterns (MPs) which are macroscopically 
different 1^ Diversity of macroscopic observable magni- 
tudes arising from random arrangements of microscopic 
constituents is the hallmark of complexity i^°d^ 

However, no consensus has yet been reached on 
whether the macroscopic limit of the Edwards- 
Andersoni^ (EA), in which only nearest-neighbor spins 
interact, (i) follows closely the SK model^^ (ii) deviates 
from SK model behavior but nevertheless shows some 
diversityji^ or (iii) fits a radically different picture, the 
droplet scenariofi^^ in which SGs with up-down sym- 
metry can only be found in one of two macroscopic spin 
configurations which are related by global spin inversion. 
Thus, according to the droplet theory, the two necessary 
ingredients (quenched disorder and frustration) for com- 
plexity would become unable to generate diversity in the 
macroscopic limit of EA systems. 

Complexity would make SGs rather exceptional among 
the many-particle systems of statistical physics. On the 
other hand, they would share this property with systems 
one finds elsewhere, such as in the life sciences?^ informa- 
tion systemSfi^ optimization problems^i^ and financci^ 
In everyone of these fields complexity and diversity are 
the rule rather than the exception. Thus, addressing 
these issues in SGs with standard methods of statistical 



physics can lead to insight into other seemingly discon- 
nected areas of research. 

The basic tool for the characterization of the SG state 
is the spin overlap q between two system statesii To de- 
fine it, let cr|^^ be the spin of system state 1 at site i, and 
similarly for 2. Then, q = N~^Yiiaf''^ af'\ i.e., q is the 
average (over all sites) spin alignment between states 1 
and 2. One usually lets states 1 and 2 be either (i) of a 
given time evolution of a given specimen at two widely 
different times?^ or (ii) of two independent (1 and 2) 
time evolutions of the same specimen. 

In macroscopic SK systems, the probability density 
(PD), p{q), averaged over all realizations of quenched dis- 
order (RQD), fulfillsii^, (i) p{q) oc f{q) -f 5{q ~ g,„) for 
g > 0, where /(q) is a smooth function of q for q < q^ 
and /(q) = for q > q^, and (ii) p(g) = p{—q) if no mag- 
netic field is applied (which we will assume throughout). 

The fact that f{q) ^ in the SK model implies the 
existence of "odd" MPs (in addition to a pair of ordinary 
MPs that one expects to observe in all magnetic systems) , 
whence complexity follows. Understandably, attempts at 
discerning between the macroscopic behaviors of the SK 
and EA models have focused on p{q). 

Specific system samples are interesting to examine. At 
least for finite-size EA systems, seed-dependent spin con- 
figurations do appear, much as in the SK modeh^, in 
thermal equilibrium. This is illustrated in Refs. [2^ and 
[23I , where plots oi pj (q) vs q are shown for two different 
J' sets of spin-spin couplings. Whereas some portions of 
Pj{q) differ drastically from sample to sample, the por- 
tions for larger values of q are alike and all of them peak 
near q = qm- We let self-overlap spikes (SOS) stand 
for spikes centered near g™. Because they are all cen- 
tered near the same position, the average of SOS (one 
for each RQD) over all RQD gives rise to the large peak 
at q ~ q„i. Spikes centered on smaller q values, which 
vary randomly with different RQD, come from spin over- 
laps between states that belong to different basins of at- 
traction. Accordingly, we refer to them as cross- overlap 
spikes (COS). 
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Interesting as it might be, statistical information on 
spike behavior has not, as far as we know, been available. 
Very little information on COS follows from the (average) 
behavior of p{q). Cross-over spike statistics would enrich 
our picture of the SG state, somewhat as the pair corre- 
lation function does for the physics of liquids. 

We aim to show how COS in pj{q), of which the lo- 
cations and shapes vary randomly over different RQD, 
can be added in a coherent fashion, in order to obtain a 
pair correlation function, g{q \ 0, Q), which is an average 
(in a sense which is defined below) over all RQD of all 
COS in the < q < Q range. We also (numerically) 
calculate the width of g{q \ 0,(5), which is a measure 
of MP thermal fluctuations, for Q = 1/2. The results 
wc obtain for low temperature (T) point to the follow- 
ing behavior: (i) g{q \ 0, 1/2) closely follows a Levy-flight 
like distribution)^ (ii) the pair correlation functions that 
follow from SOS and COS are roughly equal, that is 
giq I 0, 1/2) « g{q \ 1/2, 1), (in) the width of giq \ 0, 1/2) 
varies little, if at all, with linear system-size L (scales as 
1/L) in the EA (SK) model. Thus, different ranges of 
spin-spin interaction give rise to qualitative differences 
between the complex behavior of spin glasses. 

Models — We study the SK and EA models. In both 
of them, a (Ti = ±1 spin is located at each i-th site of 
a simple cubic lattice of N = sites. The interaction 
energy between a pair of spins at sites i and j is given 
by JijUicrj. We let Jij = ±1/\/N randomly, without 
bias, for all ij site pairs in the SK model. The transi- 
tion temperature T^g between the paramagnetic and SG 
phase is Tgg = For the EA model, Jy = imless 
ij are nearest-neighbor pairs, and we draw each nearest- 
neighbor bond Jij independently from unbiased Gaussian 
distributions of unit variance. Then, Tsg ~ 0.95i^ 

We let {u j) J stand for the average of a thermal equi- 
librium quantity u j over a number A^^ of different sets 
of random bonds {J}- 

Pair correlation function — Aiming for statistical infor- 
mation on COS at low temperature, we let 

rQ2-h{q) 

Gj{q I Qi, (32) = / dqi pj{qi)pj{qi + g), (1) 

jQl+h(~q) 

where h{q) = if g < and h{q) = q ii q > 
0. Clearly, (AQ- | q \)-^Gjiq \ Qi,Q2), where 
AQ = Q2 - Qi, is the average of pj{qi)pj{qi + q) 
over the (Qi -I- /i(— g), Q2 — h{q)) domain. We term G{q \ 
Qi,Q2) = {Gj{q I Qi,Q2))j po-if correlation function. 

The integral in Eq. ([1]) is as for the PD to be at — q 
after a two-step random walk which starts at the ori- 
gin, in which the length of both steps is identically dis- 
tributed, but they are taken in opposite directions. Note 
that Gj{q \ Qi,Q2) (i) peaks at q ~ 0^^ (h) is even 
with respect to g = 0, since pj{q) = pj{—q), and (iii) is 
somewhat broader than pj{q) (from the theory of ran- 
dom walks). 

The operation defined in Eq. ([T|) clearly displaces to 
g = any spike in pj within the (Qi, Q2) domain. Thus, 



TABLE I. The number of samples Ns and the time Ts taken 
for equilibration as well as for subsequent averaging is given 
in thousands of MC sweeps. Acceptance rates for system con- 
figuration exchanges at all T are larger than or approximately 
equal to A. AT is the temperature spacing between systems 
in the tempered MC setup. 





SK 


EA 


L 


4 6 8 


4 6 8 10 


Ts 

Ns 
AT 
A 


50 50 100 
20 10 10 
0.05 0.05 0.04 
0.7 0.5 0.4 


10 10^ 103 10^ 
30 40 30 5 
0.1 0.1 0.05 0.04 
0.7 0.4 0.5 0.45 



by appropriate choice of Qi and Q2 values, G{q \ Qi, Q2) 
enables one to make comparisons (see below) on equal 
footing of statistical information on SOS and COS. 

We can also define g{q \ Qi,Q2) = BG{q \ Qi,Q2), 
where B = 1/ Jq^Zq^ dq G{q \ Qi,Q2)- Note that 
g{q I Qi,Q2) is the (conditional) PD that 52 ~ gi = q, 
given that gi, 52 G (Qi, Q2). More specifically, assume all 
realizations of identical pairs of quenched disordered sys- 
tems are evolving independently in equilibrium. Then, 
g{q I Qi,Q2) is the PD for q = q2 — qi, in the set of 
pairs of identical systems in which (71,52 £ {QitQ2) are 
observed at infinitely far apart times. 

We define widths of two correlation functions. For a 
distribution function F{x) such that dx xF{x) = 
0, it makes sense to define a width 5x by 6xF{Q) = 
S^oa '^■^ ^(■^)- Since g{q \ Qi, Q2) is normalized, we let 

w{Qi,Q2) = l/g{0\Qi,Q2). (2) 

For short, we let w+ = w{0, 1/2), w- = w{l/2, 1). For 
T < 0.3, W- and w+ are widths for COS and SOS, re- 
spectively. Furthermore, note (i) g{q \ Qi,Q2) < 5(0 | 
Qi,Q2) implies < 1/2, (ii) that we can think of 

W- as an intrinsic width of g{q \ 0, 1/2), if w- <?C 1/2, 
and similarly for w+. We also define half- widths at half- 
maxima r_ and r+ by 2g{r_ \ 0,1/2) = g{0 \ 0,1/2) 
and 25(r+ | 1/2, 1) = ^7(0 | 1/2, 1). 

Method — All numerical results given below follow from 
parallel tempered Monte Carlo (MC) simulations 
We give all times in terms of MC sweeps. 

All pairs of systems start running from independent 
random spin configurations. Each system pair is then 
allowed to come, in time Ts, to equilibrium with each 
reservoir of a string of them at T, T + AT, T + 2AT . . ., 
before readings of q values are taken over an additional 
time span. From many such readings, the thermal equi- 
librium probability pj{q), for a given RQD, is obtained 
for each temperature. 

Relevant parameters for the simulations are in Table 

I. 

Results for the pair correlation function — Data points 
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FIG. 1. (a) Plots of g{q \ 0,1/2) vs q for EA systems at 
T = 0.2, with the values of L shown. The full line is for a fit 
to the L = 8 data points with r^/[K;(r^+ | q \'<)]. (b) Same 
as for (a) but for g{q \ 1/2, 1). (c) Same as for (a) but for the 
SK model, (d) Same as for (a) but for g{q \ 1/2, 1) for the SK 
model. Values of w and F are given in Fig. [2] for all T and L 
above, in both the EA and SK models. Values for 7 are given 
in the text. 



for the pair correlation function are shown in Fig. [T] for 
T = 0.2. The two graphs on the left hand side (right hand 
side) are for the EA (SK) model. The two top (bottom) 
graphs are for Qi,Q2 = 0, 1/2 (Qi, Q2 = 1/2, 1), which, 
since for T < 0.3, arc for COS (SOS). While neither 
COS nor SOS exhibit significant size dependence in the 
EA model, they clearly do so in the SK model. We return 
to this point below. 

All curves shown in Fig. [T] are rather pointed at the 
top. This is in contrast with the well known curves 
for p{q) in the neighborhood of g = qm for finite SK 
systems^S [but see Ref. [13 andlH for some pj{q)]. This is 
because values of qm vary over different RQD by amounts 
which, at least for L = 8 and T < 0.3, arc roughly equal 
to r+ for the SK model. Thus, averaging pj/ (q) over all 
RQD gives rise to a rounded p{q) while g{q \ 1/2, 1), be- 
ing a coherent like superposition of spikes over different 
RQD, reveals their individual shapes. 

Widths — We note that if z = xi + 2:2, and xi and X2 
are drawn from the above distribution, then, within a few 
per cent, in obvious notation, F^ ~ Y — 1 if 1.15 < 
7 < 2. Thus, g{q \ Qx.Q-i) is approximately 2/^/7 — 1 
times wider than spikes in the (Qi, Q2) domain. 

Good fits to all plots in Fig. [T] are provided by 
+ I q p)], which closely follows a Levy fiight 
distribution^, for 1 < 7 < 2. Fits to the L = 8 data 
are shown in Figs. Hlja), [Hb), [TJc) and[TKd). For both 
the EA and SK models, 7 ~ 2(1 - 1/L) for all T < 0.3. 
Values for w and F are given in all panels of Fig. [2l As in 
Fig. [U the two graphs on the left hand side (right hand 
side) are for the EA (SK) model. The two top (bottom) 
graphs are for Q\,Q2 = 0,1/2 ((5i,(52 = 1/2,1), which 




FIG. 2. (a) Log-log plots of and 2F_ vs T, for EA systems 
of 1/ = 4, 6, 8 and L — 10, as shown. Error bars for L = 10 are 
shown, but for smaller L values they are hidden by symbols, 
(b) Same as in (a) but for w+ and 2F+. (c) Same as in (a) 
but for SK systems of L = 4, 6 and L = 8, as shown, (d) 
Same as in (c) but for Lw+ and 2LF+. 




FIG. 3. Plots of (Sp^) and of {Sp^fS/L vs T for the EA and 
SK models, respectively, (5 = 1/2 and the values of L shown. 
Except for L — 10, icons cover all error bars. 



for T < 0.4 arc for COS (SOS). 

Widths W- and w+ appear in Figs. [5^ and[2jD to be size 
independent for < T < 0.3. This points to finite widths 
for COS and SOS in the L ^ 00 hmit of the EA model 
at low temperature. On the other hand, ~ 1/L 

for large L in the SK model seems consistent with the 
data points shown in Figs. [2j; and[2]i for w_ and 

Probability fluctuations over different RQD — 
Additional information follows from the (unnormalized) 
pair correlation function G{q \ Qi,Q2)- For instance, 

G{Q\Q,Q) = Q[{p'^f + {^f], (3) 

where jP and {Sp'^)^ are the averages of p and (Sp)^ over 
the <| g |< Q range, respectively. 
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Plots of {5p ) vs T for Q = 1/2 and various systems 
sizes of the EA and SK models are shown in Fig. [S] Note 

that whereas {Sp^) scales as ~ y/L in the SK model, it 
appears to be, within statistical errors, independent of L 
in the EA model. 

The amount of fluctuations over different RQD of 
Xj[Q), defined by Xj{Q) = dq pj[q), differs quali- 
tatively from {Sp^Y if Q ^^ot too small. This is because, 
{Xj{Q)) J = u;(0, Q)G(0 I 0, Q), as follows from Eq. (g]). 

Therefore, our result that w{0,Q){6p'^)^ is independent 
of L in both the EA and SK models if Q » w{0, Q) im- 
plies {Xj{Q)) J remains then bounded in both models as 
L — > oo. 

Conclusions — We have given a recipe for a coherent 
addition of self- and cross-overlap spikes (SOS and COS). 
The latter are centered on positions that vary randomly 
with RQD. In both the EA and SK models, the correla- 
tion functions for SOS and COS turn out to be approx- 
imately equal. The widths W- and r_ (both for COS) 
give a measure of the thermal fluctuations of magnetic 
patterns. They are not too different from iy+ and r+ 
(both for SOS), respectively. Neither w± nor T± vary 
much with linear system-size L in the EA model but scale 



approximately as 1/L in the SK model. Their variation 
with system size at low temperature suggests they van- 
ish in the macroscopic limit of the SK model but remain 
finite in the EA model. Finally, the mean square devia- 
tions oipj away from p{q) appear not to vary much with 
L in the EA model but scale approximately as L in the 
SK model. 

The rule we have uncovered -which relates thermal 
fluctuations of magnetic patterns as well as probability 
fluctuations to interaction range- may well be valid in 
some broader domain, beyond the SK and EA models. 
For one, preliminary (unpublished) work yields similar 
results for some spatially disordered systems which are 
geometrically frustrated. Extensions to other fields of 
complex systems easily come to mind. 
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